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The cross-plane thermal conductivity of a type II InAs/GaSb superlattice (T2SL) is measured 
from 13 K to 300 K using the 3cj method. Thermal conductivity is reduced by up to 2 orders of 
magnitude relative to the GaSb bulk substrate. The low thermal conductivity of around 1-8 W/m-K 
may serve as an advantage for thermoelectric applications at low temperatures, while presenting a 
challenge for T2SL quantum cascade lasers and high power light emitting diodes. We introduce a 
power-law approximation to model non-linearities in the thermal conductivity, resulting in increased 
or decreased peak temperature for negative or positive exponents, respectively. 

PACS numbers: 73.21. Fg,73.50.Bk,73.61.Ey,71.70.Fk,71.70.Gm 



InAs/GaSb type II broken-gap superlattices (T2SL) 
have been successfully developed for use in long 
wavelength (LWIR) and midwave (MWIR) infrared 
detectors^ - — with band engineered cut-off wavelengths 
and a reduction in dark current compared to direct gap 
bulk semiconductors. High-power LWIR and MWIR 
cascaded LEDs and lasers based on T2SL band struc- 
ture engineering^— also offer tunable infrared light emis- 
sion. Both detectors and lasers perform best at low 
temperatures^ and heat needs to be carried away from 
the devices to reduce Auger recombination and increase 
carrier lifetime for detectors, and to decrease lasing 
threshold for lasers. It is therefore important to have 
a good understanding of the thermal conductivity of this 
artificial material, to understand performance parame- 
ters at low operating temperatures. Recently, efforts by 
the authors have also identified the T2SL as a candidate 
material for low-temperature Peltier coolingii based on 
the large Seebeck coefficients that have been reported^ 
but measurements of the thermal conductivity at cryo- 
genic temperatures were lacking. In this work, we mea- 
sure the cross-plane thermal conductivity of T2SLs from 
13 K to 300 K using the 3w method ; 13 ' 14 and present a 
power-law approximation for modeling thermal conduc- 
tivity over large thermal gradients expected to occur in 
such low-thermal conducting materials. 

The two different T2SLs studied in this work were 
grown by molecular beam epitaxy on GaSb substrates. 
Because this material was intended for p-i-n detectors, it 
consists of a 0.5 /im GaSb p + (~ 10 18 cm -3 ) buffer layer, 
followed by a 0.5 /jm T2SL p + (~ 10 18 cm" 3 ) region, a 
2 fim undoped T2SL layer, a 0.5 /zm T2SL n+ (~ 10 18 
cm~ 3 ) region, and a 10 nm Si doped InAs n + capping 
layer. T2SL-1 is composed of 12 monolayers (ML) of 
InAs and 8 ML of GaSb per period. T2SL-2 is composed 
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FIG. 1: Temperature change AT vs. frequency / measured 
with the 3lo method for the T2SL + Si02 + GaSb substrate 
(filled symbols) and for Si02 + GaSb substrate (open sym- 
bols), (a) at 15 K and 300 K, and (b) at 80 K and 212 K. 
Inset: picture of the sample, with left half etched to sub- 
strate, and right half unetched. Two Au heater-thermometer 
lines are deposited on each side, with current contact pads at 
top/bottom and voltage contacts in between. 



of 19 ML of InAs and 18 ML of GaSb per period^ 

Following standard 3w sample preparation, the ther- 
mal conductivity of a sample including the layer of in- 
terest is measured relative to a reference sample without 
the layer. Thus the T2SL is wet-etched away from half 
the sample with a solution of citric acid and phosphoric 
acid plus peroxide. An insulating Si02 layer of 150 nm 
is deposited using plasma enhanced chemical vapor de- 
position (PECVD) to prevent an electrical short circuit 
through the conducting substrate. Then 200 nm thick 
gold heater-thermometer filaments are deposited atop a 
3 nm Ti adhesion layer using e-beam evaporation on both 
etched and unetched regions. The filament shown in the 
inset of Fig. (Ub) is 3.6 mm long and 30 /im wide, much 
wider than the 3 /im T2SL thickness so that the heat 
flow through the T2SL obeys the one-dimensional ther- 
mal diffusion equation* 1 ^ 

According to the 3u> method ) 13 ' 14 by sending a current 
excitation at frequency uj through the gold filament, heat- 



ing power P is induced at frequency 2u, leading to ther- 
mal variation AT at the same frequency and ultimately 
a 4-point voltage at frequency 3ui proportional to 
the temperature difference via the thermal coefficient of 
Au3 The thermal difference AT = AT sub (w) + AT f can 
have a frequency dependent substrate contribution and 
a frequency independent thin-film contribution, 



AT sub (w) 
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where K S ub is the thermal conductivity of the GaSb sub- 
strate, w and / the width and length of the Au filament 
heater, and C is a constant offset which includes the SiC>2 
layer contribution. For the T2SL film, K{ and t are the 
film thermal conductivity and thickness, respectively. 

We measured the sample in an Oxford variable temper- 
ature insert (VTI) helium gas flow cryostat from 300 K 
down to 13 K, using standard lock-in techniques. Since 
there are background 3a; voltages from the lock-in power 
source and nonlinear components in the measurement 
circuit and in the lock-in A-B input channels, we also 
measured a reference background 3cj signal with a low 
thermal coefficient resistor of equal resistance as the gold 
filament. 

By measuring the slope of the temperature difference 
AT as a function of log-frequency^, we can deduce the 
substrate thermal conductivity and thereby confirm the 
reliability of our measurements. In Fig. [21 we com- 
pare our measured GaSb substrate thermal conductivity 
(solid circles) with previously published GaSb bulk ther- 
mal conductivity (open circles) 16 ! 17 , indicating excellent 
agreement. 

The T2SL cross-plane thermal conductivity K{ is de- 
duced from the thermal difference ATf, plotted with solid 
triangles in Fig. [5J As expected from other studies of su- 
perlattice thermal conductivities, the T2SL value is re- 
duced by 2 orders of magnitude compared with the bulk 
substrate thermal conductivities for GaSb bulk (open cir- 
cles). We note that the suppression is much greater, up to 
as much as 3 1/2 orders of magnitude, when compared 
with InAs bulk (open squares). To estimate the T2SL 
thermal conductivity at arbitrary temperatures, we pro- 
vide the equation for the solid-line empirical fit to the 
measured T2SL thermal conductivity in Fig. [5] 



logK = \ogn m -B 







2 






D 


HI! 


+c 


log 















1 -i 



(3) 



with maximum thermal conductivity n m — 6.954 W/m-K 
at T m = 74 K, B = 1.0312, and C = 0.53042. 

To improve device performance, it is useful to model 
thermal distribution within the active layer of a quan- 
tum cascade laser (QCL) or LED±&±2. With poor ther- 
mal conductivity, one can expect large thermal gradients 
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FIG. 2: Cross-plane thermal conductivity for T2SL samples 
(solid triangles). Measured GaSb substrate thermal conduc- 
tivity is shown in solid circles. Published data for bulk ther- 
mal conductivity of InAs 16 and GaSb— are shown for com- 
parison in open squares and open circles respectively. The 
solid line is a fitting polynomial for the average value of T2SL 
thermal conductivity described in the text. 



across the T2SL layer, so it is useful to develop an an- 
alytical estimate of the thermal profile including large 
temperature drops. We do so below by introducing a 
power-law approximation to the temperature dependence 
of the thermal conductivity. 

The thermal profile can be determined from the power 
density generated per unit volume for the various devices 
of interest, namely LEDs and QCLs. For an infrared LED 
the power density due to Joule heating is 



P = £J 



(4) 



where the electric field £ across the T2SL emitting layer is 
assumed uniform, J is the current density and the light 
output power is neglected. For continuous-wave (GW) 
lasers, one can estimate the heat dissipated at threshold 
assuming that almost all input power is dissipated as heat 



P = (E 32 + A)J th /eL p , 



(5) 



where J t h is the threshold current density, T32 is the en- 
ergy of the optical transition, A is the energy separation 
between the ground state of the injector to the lower laser 
level of the previous active region, and L p is the length 
of one period of QCL active region and injector.—. 

Using Eqs. ((U)-©, the thermal profile in the active 
region becomes: 



P = V ■ (-kVT). 



(6) 



For a substrate-side mounted device, the primary heat 
dissipation is through the GaSb substrate, modeled as 



3 



a ID thermal diffusion problem 18 . If, for the moment, a 
constant thermal conductivity k is assumed, the solution 
to Eq. ((6]) takes the simple form~ 



. rT1 , \ 1 P 2 dT 
l kq ax 



(7) 



with boundary conditions AT(0) = and dT/dx(0) = 
dT/dx\o- 

However, we arrive at a more accurate estimate of the 
thermal profile at high powers if we include the tempera- 
ture dependence of k(T). By inspection of the T2SL log- 
log plot in Fig. [21 for large ranges of temperature above 
150 K and below 30 K the slope is roughly constant and 
k(T) can be approximated with a power-law. We thus 
empirically fit the T2SL thermal conductivity with the 
local power- law expression n(T) = ko(^) 3 , where s is 
the power-law exponent, T is the base-line temperature 
of interest, and kq is the thermal conductivity at that 
temperature. Solving Eq. (J6j> under this local power-law 
assumption, the exact solution for the temperature pro- 
file in the active T2SL region with the same boundary 
conditions becomes: 



T{x) = To 



1 + (s + 1) 



AT(x) 



T 



(8) 



where AT(x) is the same expression from Eq. ([7]) and s 
is the power-law exponent. In real devices, the change in 
temperature AT(x) will normally not exceed the absolute 
base temperature To, so this expression can be expanded 
for small AT(x) /T < 1 



T(x) ~ T 



1 AT 2 (x) 

Al (x) s 

1 ' 2 T 



(9) 



The first two terms describe the thermal distribution for 
constant baseline thermal conductivity kq, and the third 
term accounts for the local power-law assumption, pro- 
portional to the exponent s. Note that the position of 
maximum temperature does not change between Eqs. ([7]) 
and ©, nor does the thermal derivative boundary condi- 
tion at x = 0. The resulting thermal distribution is plot- 
ted in Fig. [3] under high power dissipation P such that 
the maximum temperature is about a factor of 1.5 greater 
than the substrate absolute temperature To. When the 
exponent is positive the peak temperature is reduced, for 
example s ~ 1 for T below 30 K. And when the exponent 
negative the peak temperature is increased, for example 
s 1/2 for T above 150 K. 



The local power-law exponent s in the neighborhood 
of an operating temperature To can be determined from 
the same empirical fit parameters as for Eq. ([3]): 
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FIG. 3: Calculated temperature profile in the active region 
normalized to the cold-sink substrate temperature To with 
a typical power dissipation in a T2SL LED. The power-law 
exponent s = approximates the behavior in the intermediate 
temperature range To = 50 — 100 K whereas the s — —1/2 
negative exponent is appropriate to higher temperatures To = 
150 — 300 K and s = 1 to lower temperatures To = 13 — 30 K. 



We remark that these materials could be promising as 
cryogenic thermoelectrics because below 20 K the ther- 
mal conductivity is quite low of order 1 W/m-K, and be- 
cause high Seebeck coefficients up to 2 mV/K have been 
reported for the hole band and 300 /iV/K for the electron 
band^ at 4 K. We thus conclude that careful thermal 
modeling with the knowledge of the thermal conductiv- 
ity reported here can improve device performance for a 
range of applications. 
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